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Resonances of the (Frobenius-Perron) evolution operator V for phase-space densities have recently 
attracted considerable attention, in the context of interrelations between classical and quantum 
dynamics. We determine these resonances as well as eigenvalues of V for Hamiltonian systems with 
a mixed phase space, by truncating V to finite size in a Hilbert space of phase-space functions 
and then diagonalizing. The corresponding eigenfunctions are localized on unstable manifolds of 
^~ > ' hyperbolic periodic orbits for resonances and on islands of regular motion for eigenvalues. Using 

information drawn from the eigenfunctions we reproduce the resonances found by diagonalization 
through a variant of the cycle expansion of periodic-orbit theory and as rates of correlation decay. 



■ I. INTRODUCTION 
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f""^ \ Chaotic behavior in classical Hamiltonian systems manifests itself in an infinite hierarchy of phase-space structures 
and instability of trajectories with respect to small changes in their initial condition. If resolution in phase space is 
restricted, e.g. due to finite precision of measurements, phase-space trajectories can be followed only approximately. 
Every finite-resolution initial condition actually stands for an entire ensemble of trajectories with possibly very different 
long-time behaviors. That fact suggests a probabilistic approach to the dynamics, in which a phase-space density and 
its propagator V, the so called Frobenius-Perron operator, are studied. 

Because of Liouville's theorem V can be represented by an infinite unitary matrix in a Hilbert space of phase-space 
functions. The spectrum of V thus lies on the unit circle in the complex plane. Discrete eigenvalues represent regular 
[~^. \ dynamics, while chaos entails a continuous spectrum, and along with the latter resonances of V may occur that 
characterize effectively irreversible behavior [|l]j^]. Apart from their role as decay rates, Frobenius-Perron resonances 
have been found to link classical with quantum dynamics because they could as well be identified from quantum 
systems ||[|]. Furthermore, they have been predicted to carry information on the system-specific corrections to 
universality of quantum fluctuations . 

To identify resonances, mathematically defined as singularities of the resolvent of V in some higher Riemannian 
sheet, it is necessary to analytically continue the resolvent across the continuous spectrum of V from the outside to 
the inside of the unit circle . This becomes increasingly difficult as the dynamics gains complexity. Especially 
for the large class of systems with a mixed phase space an effective approximation scheme is needed that is free of 
restrictions of previous investigations 1 10 such as hyperbolicity, one-dimensional (quasi-) phase space, or isolation 
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of the phase-space regions causing intermittency. 

In the following we present such a scheme and illustrate it for a prototypical dynamical system with a mixed phase 
space, a periodically kicked top. We represent V in a Hilbert space of phase-space functions as an infinite unitary 
matrix and subsequently truncate it to a finite A-dimensional matrix J>( N \ with resolution in phase space as the 
truncation criterion. We thus generate a nonunitary approximation to the unitary V that becomes exact as N — ► oo. 
Looking at the spectra of 'pW with increasing N, we find some eigenvalues persisting in their positions, either 
almost on the unit circle or well inside, and we focus on these "frozen" eigenvalues. While frozen (near)unimodular 
eigenvalues of turn into unimodular eigenvalues of V as N — > oo, the nonunimodular frozen eigenvalues are 

not part of the spectrum of V in that limit. We will argue that at finite but large N they rather indicate the 
positions of resonances of V . Accordingly, the eigenfunctions of are of very different nature for these two 

cases. We find eigenfunctions to unimodular eigenvalues localized and approximately constant on islands of regular 
motion in phase space that are bounded by invariant tori. In contrast, eigenfunctions to nonunimodular eigenvalues 
of V^ N ' are localized near unstable manifolds of hyperbolic periodic orbits. Just as nonunimodular eigenvalues are 
not in the spectrum of the Hilbert-space operator V, the eigenfunctions are no Hilbert-space functions in the limit 
of infinite resolution. A comparison of the eigenfunctions of P^' for different values of N and furthermore with the 
eigenfunctions of the truncated inverse operator p -1 ^) illustrates why they lie outside the Hilbert space. The strong 
localization of these (resonance-)eigenfunctions at finite N allows us to identify the relatively few groups of periodic 
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orbits associated with a particular resonance up to a certain length. From only these orbits we recover resonances via 
a variant of the so-called cycle expansion of periodic-orbit theory that works even in the case of a mixed phase space. 
To verify that frozen eigenvalues have indeed the physical interpretation of resonances, we conclude with a numerical 
experiment in which the decay of a two-point correlation function is observed. In order to find a certain resonance 
responsible for correlation decay, the initial phase-space density is chosen localized in the phase-space region where 
the resonance-eigenfunction has large amplitudes. 

A short account of the results to be discussed here has been published recently p2| . 



II. LIOUVILLE DYNAMICS ON THE SPHERE 



As a prototypical Hamiltonian system with a phase space displaying a mix of chaotic and regular behavior we 
consider a periodically kicked angular momentum vector 

( J x , J y , J z ) = (j sin 9 cos ip, j sin 9 sin tp, j cos 9) (1) 

of conserved length j, also known as the kicked top JDjfl . Such a system has one degree of freedom, and its phase 
space is a sphere. A phase-space point X = (p, q) may be characterized by an "azimutal" angle tp as the coordinate 
q and the cosine of a "polar" angle 9 as the conjugate momentum p. The dynamics is specified as a stroboscopic 
area preserving map, X' — M(X). We chose M to consist of rotations R z {(3 z ),R y {(3 y ) about the y- and z-axes and a 
"torsion", i.e. a nonlinear rotation T z (t) = R z (t cos 9) about the z-axis which changes (p by tcos9, 

M = T z {T)R z {f3 z )R y (f3y). (2) 

In the sequel we keep (3 Z — (3 y — \ fixed and vary the torsion constant t to generate a phase space which is intcgrable 
for t — and increasingly dominated by chaos as r grows. We focus on r = 2.1 as a moderately chaotic case and 
t = 10.2 as strongly chaotic (see figures |l|a,b), where elliptic islands have become so small that they are difficult to 
detect. 

In the Liouville picture the time evolution of a phase-space density p is governed by Liouville's equation, 

d tP = £p = {H,p}, (3) 

where the Liouville operator C, the Poisson bracket with the Hamiltonian H, appears as generator. For our rotations 
and torsion we separately take H = (3 y J y ,H — fl z J Zl and H = (r/2) J 2 . Denoting the respective Liouvillians by £r v , 
Cr z , and Ct z we imagine Liouville's equation for each of them integrated over a unit time span. The product of the 
resulting three propagators yields our Frobenius-Perron operator V = exp(£^) exp(£fl z ) ex.p(Cn y ). The action of V 
on the phase-space density p(q,p) can be represented in terms of the "Hamilton-picture" map M as 

p n+1 (X)=V p n (X) = e c ^e c ^e CR y p n (X) 

= J dX' S(X - M(X')) Pn (X') . (4) 

Since the map M is invertible and area preserving, det dM(X) / dX = 1, the integral simplifies to 

Vp(X)=p[M-\X)] . (5) 

While a phase-space density is usually considered as Z^-integrable it can be assumed to belong to a Hilbert space 
of L 2 -functions as well. The additional structure of the Hilbert space can be exploited to represent V by an infinite 
unitary matrix, the unitarity (Vpi, V P2) = (pi, P2) following from fla) and area preservation in the map M. The basis 
functions may be chosen as ordered with respect to phase-space resolution. This allows us to eventually truncate V 
to finite size and still retain the propagator intact for functions on large phase-space scales. 

On the unit sphere the spherical harmonics Yi m (ip,cos9) with I = 0,1,2,... and integer \m\ < / form a suitable 
basis. Phase-space resolution is characterized by the index I: using all Yi m with < / < ( max phase-space structures 
of area oc 1/Z^ax can ^ e resolved, and the number of these basis functions is N = (Z max + l) 2 . 

To check on the resolution achievable by spherical harmonics consider a Gaussian G((p, 9) — 
exp (—i(<p — ip ) 2 / A<p 2 — \(9 — 9 ) 2 /A9 2 ) with the angular widths Atp,A9 and the effective area ~ AtpA9. To 
represent it by a linear combination of spherical harmonics we need only the Yi m with |m| < too ~ 2n/Aip and 
I < 2ir/Aip + ir/A9; all other spherical harmonics have negligible scalar product with G(ip,6). Indeed, the effective 
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"wavelengths" of Yi m in ip and 9 are, respectively, 27r/|m| and ir/(l — \m\), the latter given by the number (I — |m|) 
of zeros of the Legendre function P; m in the interval < 6 < n; once these wavelengths become smaller then the 
respective widths Atp, A9 of the Gaussian the overlap of G and Y\ m is negligible. Inverting that statement we may 
say that a basis set of spherical harmonics cut off at a certain l max is suitable for representing objects on the unit 
sphere with the angular sizes ~ l maK and area <~ ^ax- 
To obtain the matrix elements 



dcos6» J dip Y* m {ip, cos 6)Y Vm , [M^(p, cos 6)] 



(6) 



for the mixed dynamics it is convenient to first determine the Frobcnius-Pcrron matrices for the two rotations 
exp(£^), exp(£p ) and the torsion exp(£T 2 ), all of which are by themselves integrable and free of chaos. Defer- 
ring the derivation to Appendix A we here just note the three matrices. The Frobenius-Perron matrix of the z- 
rotation is diagonal, 



( exp(£flj) = Su>6 mm >exp(-im(3 z ) . (7) 

V / Im.l'ra' 

For the y-rotation, (exp(£fl H ))i m ,Vm' is blockdiagonal in I with the blocks well known from quantum mechanics as 
Wigner's d-matrices, 

(expOCflj), „ ,=*«'dw03»)- (8) 

\ / Im.l'm' 

For the z-torsion the Frobenius-Perron matrix elements are given by finite sums over products of spherical Bessel 
functions ji(x) and Clebsch-Gordan coefficients, 



(exp(£ T j) =<W(- 1 )"V( 2 * + 1 )(2*' + 1 ) 

V / lm,l'm' 



l+l' 

x Yl Hf >(mr)C^V"^ m i Io- (9) 
i"=\i-i>\ 

Finally, for the mixed dynamics V is the product of the three above matrices. Due to the Kronecker deltas involved, 
each element of V is just the product of three matrix elements, one each for the rotations and the torsion. Obviously, 
V is a full matrix coupling all phase-space scales, as it is to be expected for chaotic dynamics. We furthermore remark 
that with respect to the real basis 

Via = Y t0 

yL = ^{Yim + Y l * m ), (io) 

Vim = "^/fj^ m ~ Yhn) 

the matrix V is also real , since real phase-space densities cannot become complex under time evolution. As a 
consequence, all eigenvalues of are either real or come in complex conjugate pairs. 



III. RESONANCES AND EIGENVALUES AT FINITE PHASE-SPACE RESOLUTION 

To study resonances and (discrete) eigenvalues of V we investigate the finite-size approximants P (Ar) for different 
values of N, i.e. different resolutions in phase space. In the low-resolution subspace spanned by the first N basis 
functions replaces V identically, but it completely rejects fine phase-space structures that cannot be expanded 

in terms of these N functions. 

The truncation of V to finite size breaks unitarity. Therefore, may not only have unimodular eigenvalues 

but also eigenvalues located inside the unit circle. Furthermore, the spectrum of is purely discrete. The iV- 

dependence of the eigenvalues of yields information on the spectral properties of V . Upon diagonalizing p( N "> 
and increasing N we find the "newly born" eigenvalues close to the origin while the "older" ones move about in the 
complex plane. "Very old" ones eventually settle for good. If the classical dynamics is integrable (r = or (3 y = 0), 
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the asymptotic large-iV loci are back to the unit circle, where the full V has its discrete eigenvalues. The situation 
is different for the mixed phase space: while some eigenvalues of V^ N ' "freeze" with unit moduli, others come to 
rest inside the unit circle as N — > oo. These eigenvalues reflect resonances of V in a higher Riemannian sheet of the 
complex plane. We see this phenomenon as analogous to the spectral concentration in perturbation theory ]l4|,[r|jL5|| : 
a perturbation series for an operator with continuous spectrum does not converge but produces, with increasing order, 
a sequence of points concentrated in the neighborhood of the respective resonance. 

An additional intuitive argument for the persistence of non-unimodular eigenvalues as N — > oo for nonintegrable 
dynamics is the following: in contrast to regular motion, chaos comes with a hierarchy of phase-space structures which 
extends without end to ever finer scales. Probability that is propagated from large to fine scales mostly does not 
return. This effectively dissipative character of the conservative dynamics is accounted for as true dissipation by the 
truncated propagator and its nonunimodular eigenvalues, but the difference between effective and true dissipation is 
irrelevant for low resolution. Moreover, a frozen nonunimodular eigenvalue is to be understood as a scale-independent 
dissipation rate, in tune with the expected selfsimilarity in phase space. 

The freezing of non-unimodular eigenvalues is illustrated for two kicked tops, one with weakly chaotic (r = 2.1) 
and one with strongly chaotic character (r = 10.2). Figures [l]a and |l|b show the respective phase-space portraits. 
While for r = 2.1 large islands of regular motion still exist, for r = 10.2 regular islands have become so small that 
they are difficult to detect and impossible to resolve with l max = 70. Figures ||a and ||b show grey-shade coded 
histograms in the complex plane for the number of eigenvalues with moduli larger than 0.2 of all matrices with 
( max = 40, 41, 42, . . . 70. Eigenvalues with moduli smaller than 0.2 are rejected because they have not settled yet 
and their density near the origin is so high that they would spoil the histogram. Large amplitudes (black) in the 
histograms indicate the positions of frozen eigenvalues. While for r = 2.1 a large number of frozen eigenvalues is found 
close to the unit circle, in the case r = 10.2, where no regular islands can be resolved, the only unimodular eigenvalue 
is the one at unity which pertains to the stationary uniform eigenfunction Yqq; all other frozen eigenvalues lie well 
inside the unit circle and indicate resonance positions. Table [j] lists some of these frozen nonunimodular eigenvalues 
at resolutions i max = 30, 40, 50 and 60. 

At first sight the two cases presented above may appear to correspond to quite different types of dynamics, but 
actually this is only a matter of resolution. As resolution is further increased in the strongly chaotic case (t = 10.2), 
elliptic islands will be resolved at some stage. Correspondingly, eigenvalues oiV^ N ^ will emerge which freeze on the 
unit circle. Also new nonunimodular frozen eigenvalues may appear, possibly even ones larger in modulus than the 
frozen eigenvalues found at lower resolution. 



IV. EIGENFUNCTIONS AND RESONANCE- EIGENFUNCTIONS AT FINITE PHASE-SPACE 

RESOLUTION 

An investigation of the eigenfunctions of reveals more about the nature of the frozen eigenvalues. Eigenvalues 
freezing with unit moduli have eigenfunctions located on elliptic islands of regular motion surrounding elliptic periodic 
orbits in phase space. Such islands are bounded by invariant tori which form impenetrable barriers in phase space. 
Thus if Xp,n is a characteristic function constant on an island around the nth point of an p-periodic orbit and zero 
outside with x P ,n+i = ^Xp.n and x P ,i = 'PXp^ tnen 

p 

fp = ^2Xp,n (11) 

n=l 

is an eigenfunction of V with eigenvalue unity. Moreover, linear combinations of such functions are eigenfunctions as 
well. Figures ^a,b show two eigenfunctions of with eigenvalues 0.999976 and 0.999974 almost at unity for the 
kicked top with r = 2.1 and l max = 60. Amplitudes of the functions are large in dark shaded areas of phase space. 
The elliptic islands supporting the eigenfunctions can easily be identified in the phase-space portrait in figure [l]a. 

For p > 1 we furthermore expect and find the pth roots of unity as eigenvalues as well with eigenfunctions that 
have constant moduli on elliptic islands and are invariant under V p , 

p 

fp, m = ^ 2nmn/p X P ,n, withm = l,...,p. (12) 

71=1 

Superpositions of two such f Vtm , f p > m > are also eigenfunctions if m/p — m'/p', i.e. if all periods p' involved are 
multiples of p; the possible phases are dictated by the shortest period p. 
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Once freezing has been observed for eigenvalues with moduli smaller than unity their eigenfunctions have approached 
their final shape on the resolved phase-space scales. In contrast to the eigenfunctions with unimodular eigenvalues 
these eigenfunctions are localized around unstable manifolds of hyperbolic periodic orbits, ones with low periods at 
first since these are easiest to resolve. [Typically, for every period several neighbouring orbits of approximately equal 
instability can be identified in an eigenfunction, their number being the larger the more extended the eigenfunction.] 
Figure |^ shows moduli of eigenfunctions for some of the nonunimodular eigenvalues listed in table || at resolution 
( max = 60. Particularly in figures |^b,c the localization on orbits of length six and four respectively is obvious. 
With growing Z max more orbits of higher periods appear in the eigenfunction. Their supports gain complexity, in 
correspondence with the infinitely convoluted shape of the unstable manifolds. Just as for the eigenvalues there is no 
strict convergence of the eigenfunctions. Since no finite approximation J>( N ) accounts for arbitrarily fine structures 
one encounters the aforementioned loss of probability from resolved to unresolved scales. Not even in the limit N — > oo 
can the unitarity ofV be restored: rather, in tune with a continuous spectrum of V, the eigenfunctions tend to singular 
objects outside the Hilbert space. This becomes even more obvious in a comparison of eigenfunctions of V^ N ^ with 
eigenfunctions of the truncated inverse propagator 'P -1 '^. Both matrices have the same frozen eigenvalues, but only 
for the unimodular ones eigenfunctions of and V^ 1 ^ to the same eigenvalue have (roughly) the same support, 
as it is expected for unitary dynamics: the elliptic islands that support the eigenfunctions arc invariant under time 
reversal. For nonunimodular frozen eigenvalues eigenfunctions of and r p~ 1 ( N '> have distinct supports. While 

unstable periodic orbits are traversed in the opposite direction but remain the same under time reversal, stable and 
unstable manifolds are interchanged. Thus the eigenfunctions of P" 1 ^) are localized near the same periodic orbits 
but around the stable manifolds of the forward-dynamics. In particular the eigenfunctions do not converge to the 
same limit as N — ► oo. Needless to say this again reflects the effective irreversibility of the chaotic Hamiltonian 
dynamics. Figures || and ^ illustrate the aforesaid. The eigenfunction in figure [5]a at resolution l max = 30 corresponds 
to the same frozen eigenvalue as the eigenfunction in figure at resolution Z max = 60. While coarse structures are 
the same in both eigenfunctions, more fine structures appear at the higher resolution. Figure ||c shows, again for the 
same eigenvalue, the eigenfunction of ■p~ 1 ' Ar - ) at Z max = 60. In contrast to figure [sJd localization is now on the stable 
manifolds of the same periodic orbits of the foreward dynamics. These periodic orbits as well as their unstable and 
stable manifolds are shown in figure ^. 

Just like the unimodular eigenvalues, frozen eigenvalues inside the unit circle come in p-families, i.e. with phases 
corresponding to the pth roots of unity, determined by the length p of the shortest periodic orbit present in an 
eigenfunction. Now the above argument must be modified in the following way: assume an eigenfunction / is mostly 
concentrated around a shortest unstable orbit with period p as well as a longer one with period p 1 . Denote by ipp,n 
again a "characteristic function", now constant near the n-th point of the period— p orbit, n — l...p, and zero 
elsewhere. The truncated Frobenius-Perron operator maps into V^ N ^ p ^ n = r p tl>p,n+l with the real positive 
factor r p smaller than unity accounting for losses, in particular to unresolved scales. As previously, independent linear 
combinations of the ip p , n can be formed as 

p 

f P , m = J2 ei ^ mn/P ^ n withm = l...p. (13) 

n=l 

Now consider a sum of two such functions, g = f Ptm + f p ', m >, and apply pW, For g to qualify as an approximate 
eigenfunction we must obviously have r p « ry, i.e. the two periodic orbits must be similarly unstable, and m/p — 
ml jp 1 . But then indeed V^ N 'g ~ r p e l27rm / p g and [pwpj ss r p g, and the phase is dictated by the shortest orbit. Since 
orbits of low period are most likely to be resolved first, usually the first frozen eigenvalues found from diagonalizing 
pW have phases corresponding to small numbers p. Since nonunimodular eigenvalues with phases according to 
p = 1,2,4 have already been presented in |jT^| , in this paper we intentionally chose r = 10.2, a case for which 
eigenvalues with other phases happen to be more easily detected. In figure and table | complex eigenvalues with 
phases near the sixth roots of unity are clearly visible. 

V. PERIODIC-ORBIT EXPANSION OF THE SPECTRAL DETERMINANT 

Since the unstable periodic orbits that are linked to a nonunimodular eigenvalue can be identified quite easily from 
the eigenfunction of pW, one is tempted to adopt a cycle expansion to calculate decay rates from periodic orbits. 

For purely hyperbolic systems, cycle expansions of the spectral determinant, i.e. the characteristic polynomial of 
the Frobenius-Perron operator 
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d(z) = det(l - zV) = exptrln(l - zV) 



(14) 



are known to allow for the calculation of resonances with high accuracy |lq , |17| . The spectral determinant is expressed 
in terms of the traces of the Frobenius-Perron operator tr V n as 



d(ar)= JJexpj-— trP") (15) 



and subsequently expanded as a finite polynomial up to some order n max . Only the first n max traces are required for 
the calculation of this polynomial, and only periodic orbits with primitive periods up to n max contribute to it. 
The traces 



tr V" = J dX S(X - M n (X)) (16) 
are calculated as sums over hyperbolic periodic orbits of length n as 

where the matrix S„ = dM n (X) / dX is the linearized version of the map M n evaluated at any of the points of a 
period— n orbit. If n is the r-fold of a primitive period p, the matrix 3 n can as well be written as S„ = 3£, and the 
spectral determinant as a product of contributions from all primitive orbits {£;}, 



d(z) = J|exp 
{k} 



E 1 



^rdet(l-^J 



(18) 



The zeros z% of the polynomial which are insensitive against an increase of n max are inverses of resonances. 

For the ordinary cycle expansion of a spectral determinant just sketched to converge, all periodic orbits must be 
hyperbolic and sufficiently unstable |j],fl6|,jl7| . These conditions are not met for a mixed phase space, not only because 
of the presence of elliptic orbits but also because of very weakly unstable orbits with high periods. We can circumvent 
such limitations by considering only one ergodic region in phase space at a time, bar contributions from elliptic orbits, 
and impose a stability bound by including only the relatively few, about equally unstable hyperbolic orbits 
identified in an eigenfunction of pW. We surmize that the spectral determinant factors as 

oo 

d(z) = l[d i (z), (19) 

i=\ 

with one factor di for the family {k}i of periodic orbits showing up in the zth eigenfunction. Each such factor di(z) is 
then calculated separately with the above expressions. Obviously this variant of the cycle expansion reproduces the 
phases of the resonances exactly since these are again directly determined by lengths of orbits: if p is the smallest 
period used in di(z), the polynomial can as well be written as a polynomial in z p thus allowing the zeros to have the 
phases of the pth roots of unity. Needless to say, the notorious proliferation of periodic orbits with growing length 
causes more and more difficulties for the cycle expansion for p-families of resonances with increasing p. 

We first turn to the top with r = 10. Here, unstable orbits with periods 1,2, and 4 and the corresponding p-families 
of resonances are easily detected. In Table || frozen eigenvalues found at the resolution Z max = 60 are listed in the first 
column. The columns to the right show the corresponding results of the cycle expansion up to order n max = 1,2,4. 
The total number of orbits used in each case is given in curly brackets. For instance, our cycle expansion of the 
resonance 0.81 involves three orbits of lengths 1, 2 and 4. The resonance —0.7510 as well as its phase-related partners 
—0.0079 ± i0.7517 can be approximately reproduced using four orbits of length 4. For the real positive eigenvalue 
0.7470 (not listed in table ||) completing this quartet of eigenvalues the eigenfunction is also localized in other regions 
of phase space and thus further orbits with periods different from p — 4 would have to be included in the cycle 
expansion. The first repetition of the single period-2 orbit contributing to the resonance 0.6597 in table || gives an 
almost diverging contribution to the spectral determinant, thus hindering its expansion to a higher order. 

The top with r = 10.2 offers itself for a study of complex frozen eigenvalues with phases near the fourth and 
sixth roots of unity. Even though these two groups of eigenvalues are almost equal in modulus, the corresponding 
eigenfunctions are localized in different regions of phase space, and therefore the two examples turn out independent 
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(see figures @3,c). The reason why we leave aside the eigenvalues on the real axis is that their eigenfunctions show 
localization on many orbits of different, higher periods which cannot all be identified (see figures ||a,d). For the 
resonances —0.0058 ± i0.7080 the spectral determinant of order n max = 4, calculated from six periodic orbits of length 
4, yields ±i0.8150. In the next higher order n max = 8 another 9 orbits of primitive length 8 must be taken into 
account with the resulting approximants ±i0.7678. For the resonances 0.3550 ± i0.6199 and —0.3388 ± i0.6243 with 
phases near the sixth roots of unity 12 orbits of primitive length 6 contribute to the spectral determinant in lowest 
order. These orbits are shown in figure |^a. In lowest order the resulting six approximants of modulus 0.8599 still 
differ by about 20% from the respective moduli 0.7144 and 0.7103. These results are compiled in Table [II. 



VI. RESONANCES VERSUS CORRELATION DECAY 



To check on the physical meaning of frozen eigenvalues with moduli smaller than unity as Frobenius-Perron reso- 
nances we propose to compare their moduli with rates of correlation decay. In a numerical experiment we investigate 
the decay of the correlator 

_ (p(n)p(O)) - (p(oo)p(O)) 
C{U) - Wf) <p(oo)p(0)) ' (20) 

The initial densities p(0) are chosen as characteristic functions on a grid of 2000 x 2000 cells in phase space. Each 
cell is represented by its central point (q,p) c , and its temporal successor after one time step is the cell containing 
M((q,p) c ). After n iterations we obtain p(n). From the phase-space average (p(n)p(O)) we subtract the correlations 
remaining at n = oo due to the compactness of the phase space, and subsequently normalize. 

Depending on the choice of p(0) different long-time decays are observable. We choose p(0) as covering the regions 
where the hyperbolic orbits relevant for a given resonance are situated. The long-time decay turns out insensitive 
to the exact choice of p(0), provided p(0) does not extend to phase-space regions supporting different eigenfunctions 
with frozen eigenvalues of larger moduli. Since p(0) is real and positive, and positivity is preserved under the time 
evolution, we find positive real decay factors C(n+ 1)/C(n). Indeed, after two or three time steps these decay factors 
are in rather good agrement with the modulus of the corresponding nonunimodular eigenvalue as well as with the 
result of our cycle expansion. Tables [0] and |l| list the decay factors we obtained from a numerical fit of C(n) from 
n = 2 to n = 18, together with the moduli of the respective resonances. Figure fj] shows the decay of C(n) (dots), with 
p(0) chosen according to the eigenfunction shown in figure [|b. The numerical fit (line) yields a decay factor 0.7706, 
to be compared with 0.7103 as the modulus of the eigenvalue. 



VII. CONCLUSIONS 



In conclusion, we have presented a method to determine Frobenius-Perron resonances and eigenvalues which is 
applicable to systems with a mixed phase space. Resonances as well as eigenvalues are identified as frozen eigenvalues 
of a truncated propagator matrix in a Hilbert space of phase-space functions. The corresponding eigenfunctions 
of pW) are strongly localized on the associated phase-space structures, elliptic islands for eigenvalues and unstable 
manifolds of hyperbolic periodic orbits for resonances. The expected characteristics of the resonance eigenfunctions, 
an infinite hierarchy of structures and the distinct supports of foreward and backward time evolution, are clearly 
visible. The resonances obtained can be reproduced by a variant of the cycle expansion in which only the orbits 
that are identified in the eigenfunctions are taken into account. Finally, the correlation decay within a phase-space 
region that an eigenfunction is localized in is well predicted by the corresponding resonance. The cycle expansion as 
well as the observation of correlation decay may be used to check the accuracy of the resonances obtained as frozen 
eigenvalues. An accurate knowledge of the relevant resonances of mixed dynamics at a certain phase-space resolution 
is a desirable piece of information for many further investigations such as the role of classical resonances in quantum 
dynamics. 

We are grateful to Shmuel Fishman for discussions initiating as well as accompanying this work. Support by the 
Sonderforschungsbereich "Unordnung und grofie Fluktuationen" and the Minerva foundation is thankfully acknowl- 
edged. 
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VIII. APPENDIX A: MATRIX ELEMENTS OF THE FROBENIUS-PERRON OPERATOR 



We here derive the matrix elements of the Frobenius-Perron operator P. It is convenient to use the representation 
of Yi m in terms of associated Legendre Polynomials P[ a (x) Jl9[ . 



Yi m (q,p) = (-If 



21 + 1 (l-m)\ 



An (l + m)\ 



(21) 



For rotations R y , R z and torsion T z on the unit sphere the explicit forms for M read 



Ry(Py) 



arg — p 2 (cos {3 y cos q + i sin q) +p sin {3 y ^J 
p cos P y — y/l — p 2 sin P y cos q 



'q + Pz 
P 

■ q + rp' 
P 

The above immediately implies the intuitive result for the z-rotation 



(22) 



Y lm (R7 1 (q,p))=e- im ^Y lm ( q ,p)- 



(23) 



Due to the orthonormality of the spherical harmonics we then find the Frobenius-Perron operator represented by the 
diagonal matrix 



exp(£ fi J; imii 



hm 2 — "hh "mim2 



-im 2 /3z 



(24) 



Obviously and intuitively, exp(£# z ) is identical to the quantum mechanical (J% = 1) rotation matrix 
(ii7Tii| exp(— i/3js-L z ) 1^2*712) involving the L 2 -componcnt of the angular momentum operator L = (L x , L y , L z ) T and 
the eigenbasis Y/ m = \lm) with L 2 \lm) = 1(1 + l)\lm) and L z \lm) — m\lm). Since this formal identity holds for the 
z-rotation it must also hold for the y-rotation. Thus for exp(£fl H ) we must take the representation of the quantum 
operator exp(— if3 y L y ), well known as Wigner's d-matrix, 



exp(£fljjimi, 



(hm 1 \exp(-ip y Ly)\l 2 m2) = 5 hh d^ im2 (/3 y ), 



(25) 



which can be obtained from recursion relations p8[ . This matrix is still diagonal in l\,l2- In fact, the Frobenius-Perron 
matrix for any rotation is diagonal in that pair of indices. An immeditate implication is that pure rotations do not 
couple the subspace of Yi m with I < ^ max to its complement. For the integrable dynamics of rotation no probability 
is dissipated from the subspace and the spectrum of V^' lies exactly on the unit circle. 
For the z-torsion we make the replacement 



r /m (T : T 1 ( 9 ,p)) = e- i ™f Y lm (q,p) 
in the phase-space integral (0) and then use the expansion p5| 



-ik cos 6 



(-i)VM2i + l) 31 (k) Y l0 cos 6) , 
1=0 

which involves the spherical Bessel function ji(k). The parity Y* m = (—l) mi Y^- mi helps us on to 

OO r. 

exp(£ T Ji imi ,i 2m2 = (-l) mi V(-i)V47r(2/ + l)jj(mir) / dq dp Y lQ Y hi _ mi Y hm2 . 



(26) 



(27) 



(28) 



1=0 



We here meet Gaunt's integral [Ell over three spherical harmonics which can be written in terms of two Clebsch-Gordan 
coefficients, 



dq dp F z *„ Y hmi Y hm2 = 



(2h + l)(2i a + 1) 



4tt(2/ + 1) 



°000 °mim 2 m • 



(29) 
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Finally, using the following properties of Clebsch-Gordan coefficients, 



a 



h h I 



if l<£ {(h + l2),(h + h)-l,...,\h-l2\}, 



hhl =0 if m 1 + m 2 ^m. 



(30) 



hhi 

000 



if h + l 2 + 1 odd 



the infinite sum in (28) can be reduced to a finite one, 



h+h 



exp(£ T J iimi)i2m2 = S mim2 (-l) m W(2h + lK2l 2 + l) H)' C ooo l C -l\Lo ■ (31) 

i=\h-h\ 



It is worth noting that due to the properties (|30|) the matrix elements (J3lJ) are real if h+h is even or purely imaginary 
if h + h is odd. Furthermore they are symmetric with respect to h and h- 



IX. APPENDIX B: RECURSION RELATION FOR THE TORSION MATRIX ELEMENTS 



A recursion relation similar to the recursion relation used to calculate Wigner's d-matrix can be derived for the 
torsion matrix elements. For simplicity we consider the integral 



t z {l', l,m) = J du P^iu) e 1Tmu P/ n (w) 

which is related to the torsion matrix elements by 

~(2l' + l)(V -m)! (21 + !)(/- m)\ 



exp(£ T J 



I'mlr, 



(Z' + to)! (l + m)\ 

The associated Legendre polynomials obey the differential equation 

d 2 . d 



(!- u JTS -2u— + 1(1 + 1) - 



du 



1-U A 



tg(l',l,m) . 



pr(u) = i). 



This suggests to draw a recursion relation from 



du [DP^(u)}Pr(u)e^ 



In the case m ^ integration by parts leads to 



dupy e irmu [(2itto) ((i - u2 )^ p r - u P D + 



Now with the help of the identity 



(l'(V + 1) - 1(1 + !)-(!- u 2 ) r 2 m 2 ) PP 1 ] = 



(1 - u 2 )— Pr = -luPT + + m)PlH 1 
du 



the derivative in (|3^) can be eliminated and we arrive at the five-step recursion relation 

1 (2/ + 1)(2/ + 3) 



t z (l',l + 2,m) 



T 2 m 2 (I - m + 1) (I - m + 2) 
-(l'(l' + 1) -1(1+ 1) -2t 2 to 2 



2irTO - + 1 ^ - ™ + t z (l', l + l,m) 
l(l + l) + m 2 -l 



41(1 + 1) -3 



t z (l', I, to) 



-2irm 



1(1 + m) 
~2f 



2 2 + Tn)(l + TO — 1) 



7 fill \ 2 2v ' /\ ' 7 fill c\ \ 

— t z [l , t — 1, to) — r m — r-; 77— r — t z (t , t — 2, to) 

1 zV ' ' y (2/ + 1)(2Z - 1) v ' ' ; 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 
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For constant indices I' and m ^ this links matrix elements with five successive values of the index I. Just like 
the torsion matrix elements for these indices the terms in the recursion relation have alternately real and imaginary 
values. Even though the recursion relation involves five matrix elements, only two integrals have to be done initially. 
These are the nonvanishing integrals t z (l', \m\,m) and t z (l', \m\ + l,m) for the smallest values of I. Together with 
t z (l', \m\ — l,m) = and t z (l', \m\ — 2, m) = they can be used to start the recursion. 
For m = the recursion formula does not apply, but the integral 

**( r '*'°) = 2zTT <y "' (39) 

is trivial in this case. 
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/ — 30 

jiinx. — 


Tn q>x — 40 


^inax — 50 


/ YY1 ax 6 


0.7700 


0.7688 


0.7523 


0.7696 


0.3075 


0.3429 


0.3523 


0.3550 


±i 0.5740 


±i 0.6140 


±i 0.6211 


±i 0.6199 


-0.3170 


-0.3348 


-0.3444 


-0.3388 


±i 0.6003 


±i 0.6272 


±i 0.6283 


±i 0.6243 


-0.0042 


-0.0002 


-0.0100 


-0.0058 


±i 0.7161 


±i 0.7133 


±i 0.6930 


±i 0.7080 


-0.7025 


-0.7228 


-0.7155 


-0.7165 


0.6544 


0.6230 


0.6495 


0.6480 




-0.5619 


-0.5753 


-0.5667 



TABLE I. Frozen nonunimodular eigenvalues for r = 10.2 at Z max = 30, 40, 50 and 60 that can be identified from the 
histogram shown in figure 0b. 



^max — 60 Hmax — 1 ^max 2 


"max = 4 Decay of C(n) 


0.8103 0.2185{1} 0.7070{2} 


0.7664{3} 0.8005 


-0.7510 


-0.7483{4} 0.7697 


-0.0079 
±i 0.7517 


-0.7483{4} 0.7697 


0.6597 ... 0.4969{1} 


0.6783 



TABLE II. Left column: resonances obtained from the truncated propagator for r = 10.0 and Z max = 60. Middle columns: 
corresponding results from cycle expansion up to order n max = 1, 2, 4. The total number of primitive orbits employed is given 
in curly brackets. The eigenvalues in the second and third row belong to the same quartet of phase-related partners. Right 
column: Associated decay factors by which C(n) decreases over one time step, obtained from numerical fit. 



^max — 60 ^max — 4 


^max — 6 


n maj = 8 Decay of C (n) 






±i 0.7678{15} 0.7694 


0.3550 
±i 0.6199 


0.4299 
±i 0.7447 1 ' 


0.7706 


-0.3388 
±i 0.6243 


-0.4299 
±i 0.7447 1 ; 


0.7706 



TABLE III. Left column: resonances obtained from the truncated propagator for r = 10.2 and Z max = 60. Middle columns: 
corresponding results from cycle expansion up to order n max = 4, 6, 8. The total number of primitive orbits employed is given 
in curly brackets. The eigenvalues in the second and third row belong to the same sextet of phase-related partners. Right 
column: Associated decay factors by which C(n) decreases over one time step, obtained from a numerical fit. 
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-7t q TT 

FIG. 3. Eigenvalues 0.999976 (a) and 0.999974 (b) almost at unity of V (N) with Z max = 60 in the weakly chaotic case r = 2.1 
have eigenfunctions that are localized on elliptic islands. Dark-shaded regions in phase space indicate large amplitudes of the 
eigenfunctions. 



14 





it q k 

FIG. 4. For frozen eigenvalues 0.7696 (a), -0.3388±i0.6243 (b), -0.0058±i0.7080 (c), 0.6480 (d) of V (N) (r = 10.2, Z max = 60) 
inside the unit circle eigenfunctions have large amplitudes (dark-shaded areas) around unstable manifolds of hyperbolic periodic 
orbits in phase space. 
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(a) 



(b) 



1 TTT 



I 



m 



-n q n 

FIG. 5. As phase-space resolution increases from l n 



= 30 (a) to Zmax = 60 (b), the eigenfunction of T>^ N > for the frozen 
nonunimodular eigenvalue —0.3388 ± i0.6243 gains new structures on finer scales (r = 10.2). The corresponding eigenfunction 
(c) of J>~ l ( N } (resolution / max = 60) is localized at the same periodic orbit as the eigenfunction (b) but with stable and unstable 
manifolds interchanged. Therefore, see also figure ^j]a for the periodic orbits, figure ^|b for the unstable and figure |^c for the 
stable manifolds; these manifolds are clearly dominating the present eigenfunction. 



1G 



M 




-ir q 



FIG. 6. (a): The 12 orbits of primitive length 6 that can be identified for the resonances 0.3550±i0.6199 and — 0.3388±i0.6243 
with phases near the sixth roots of unity from the eigenfunction shown in figures ^b (also ||b). They contribute to the expansion 
of di(z) in lowest order. The unstable manifolds of these orbits are shown in (b), the stable manifolds in (c). 



17 



I i 




5 (0 15 II 

FIG. 7. Decay of C(n) (dots) with p(0) chosen according to the eigenfunction shown in figure ^> (also |E|b). The numerical 
fit (line) yields a decay factor 0.7706. 
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